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Abstract 

We study the Wigner functions of the nucleon which provide multidimensional images of the 
quark distributions in phase space. These functions can be obtained through a Fourier transform 
in the transverse space of the generalized transverse-momentum dependent parton distributions. 
They depend on both the transverse position and the three-momentum of the quark relative to 
the nucleon, and therefore combine in a single picture all the information contained in the gen- 
eralized parton distributions and the transverse-momentum dependent parton distributions. We 
focus the discussion on the distributions of unpolarized/longitudinally polarized quark in an unpo- 
larized/longitudinally polarized nucleon. In this way, we can study the role of the orbital angular 
momentum of the quark in shaping the nucleon and its correlations with the quark and nucleon 
polarizations. The quark orbital angular momentum is also calculated from its phase-space average 
weighted with the Wigner distribution of unpolarized quarks in a longitudinally polarized nucleon. 
The corresponding results obtained within different light-cone quark models are compared with 
alternative definitions of the quark orbital angular momentum, as given in terms of generalized 
parton distributions and transverse-momentum dependent parton distributions. 
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I. INTRODUCTION 



One of the most challenging tasks for unravelling the partonic structure of hadrons is 
mapping the distribution of momentum and spin of the proton onto its constituents. To this 



aim, generalized parton distributions (GPDs) [l|-|6| and transverse-momentum dependent 



parton distributions (TMDs) [7H12| have proven to be among the most useful tools. GPDs 



provide a new method of spatial imaging of the nucleon |13Ml7l|. through the definition 
of impact-parameter dependent densities (IPDs) which reveal the correlations between the 
quark distributions in transverse-coordinate (or impact-parameter) space and longitudinal 
momentum for different quark and target polarizations. On the other hand, TMDs contain 
novel and direct three-dimensional information about the strength of different spin-spin and 
spin-orbit correlations in the momentum space The ultimate understanding of the 

partonic structure of the nucleon can be gained by means of joint position-and-momentum 
(or phase-space) distributions such as the Wigner distributions. These distributions contain 
the most general one-body information of partons, corresponding to the full one-body den- 
sity matrix in both momentum and position space, and reduce in certain limits to TMDs 
and GPDs. Because of the uncertainty principle which prevents to know simultaneously the 
position and momentum of a quantum-mechanical system, the phase-space distributions do 
not have a density interpretation. Only in the classical limit they become positive definite. 
Nonetheless, the physics of a phase-space distribution is very rich and one can try to select 
certain situations where a semi-classical interpretation is still possible. Wigner distributions 
have already been applied in many physics areas like heavy ion collisions, quantum molec- 



ular dynamics, signa l analysis, quantum information, optics, image processing, non 
dynamics, . 



22|-|24f . and can even be measured directly in some experiments 



25 



-28 



mear 



The concept of Wigner distributions in QCD for quarks and gluons was first explored 



in Refs. 



29 



301 ] . introducing the definition of a Wigner operator whose matrix elements 



in the nucleon states were interpreted as distributions of the partons in a six-dimensional 
space (three position and three momentum coordinates). The link with GPDs was exploited 
to obtain three-dimensional spatial images of the proton which were interpreted as charge 
distributions of the quarks for fixed values of the Feynman variable x. This interpretation 
relies however on a nonrelativistic approximation. 

Wigner distributions have a direct connection with the generalized parton correlation 
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functions (GPCFs) which were recently introduced in Ref. 311 ] . The GPCFs are the dis- 
tributions that parametrize the fully unintegrated, off-diagonal quark-quark correlator for a 
hadron. In the case of the nucleon and after integration over the light-cone energy of the 
quark, one finds the so-called generalized transverse-momentum dependent parton distri- 
butions (GTMDs). At leading- twist there are 16 GTMDs which depend on the light-cone 
three-momentum of the quark and, in addition, on the momentum transfer to the nucleon 
A M . After two-dimensional Fourier transform from Aj_ to the impact-parameter space co- 
ordinates b±, in a frame without momentum transfer along the light-cone direction, one 
obtains the Wigner distributions which are completely consistent with special relativity. 

The purpose of this paper is to investigate the phenomenology of the quark Wigner distri- 
butions. As a matter of fact, since it is not known how to access these distributions directly 
from experiments, phenomenological models are very powerful in this context. Collecting 
the information that one can learn from quark models which were built up on the basis 
of available experimental information on GPDs and TMDs, one can hope to reconstruct 
a faithful description of the physics of the Wigner distributions. To this aim we will rely 
on models for the light-cone wave functions (LCWFs) which have already been used for 
the description of the generalized parton distributions (GPDs) jf], [32 -34], the transverse- 



momentum dependen 



of the nucleon 



34 



38 



; parton distributions (TMDs) 2l|, |34J437|| and electroweak properties 



441]. 



The plan of the manuscript is as follows. In sec. Ill Al we review the definition of the 
Wigner distributions obtained by Fourier transform of the GTMDs to the impact-parameter 
space. Although the Wigner distributions cannot have a strict probabilistic interpretation, 
they reduce to genuine probability distributions after integration over position and/or mo- 
mentum variables. As discussed in sec. Ill B\ one can obtain four types of three dimensional 
densities: in addition to momentum distributions given by the TMDs and to IPDs related 
to the GPDs, there are two new distributions mapping the nucleon as functions of one 
transverse-space coordinate and one transverse-momentum component which are not conju- 
gated and therefore not constrained by the Heisenberg uncertainty principle. Whereas the 
GTMDs are in general complex-valued functions, the two-dimensional Fourier transforms 
of the GTMDs are always real-valued functions, in accordance with their interpretation as 
phase-space distributions. These 16 functions can be disentangled by selecting different 
configurations, along three orthogonal directions, of the nucleon and quark polarizations. 
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In order to simplify the discussion, in sec. Ill CI we focus on the cases without transverse 
polarizations. In sec. Ill Dl we discuss and compare different definitions of the quark orbital 
angular momentum, as obtained from GPDs, TMDs and Wigner distributions. In particular, 
by treating the Wigner functions as if they were classical distributions, we can obtain the 
expectation value of the orbital angular momentum operator from its phase-space average 
weighted with the Wigner distribution of unpolarized quarks in a longitudinally polarized 
nucleon. In sec. II III we explicitly calculate the Wigner distributions in two light-cone quark 
models, showing the results for the first x moments in the four-dimensional phase space (two 
transverse position and two transverse momentum coordinates). In particular, we discuss 
specific situations where the density matrices have a quasi-probabilistic interpretation, giv- 
ing a semi-classical picture for multidimensional images of the nucleon. In section [IV] we 
draw our conclusions. 



II. WIGNER DISTRIBUTIONS 



A. Wigner Operators and Wigner Distributions 



Wigner distributions in QCD were first explored in Refs. 29|, |30|. Neglecting relativistic 
effects, the authors used the standard three-dimensional Fourier transform in the Breit frame 
and introduced six-dimensional Wigner distributions (three position and three momentum 
coordinates). We propose to study instead five-dimensional Wigner distributions (two posi- 
tion and three momentum coordinates) as seen from the infinite momentum frame (IMF). 



The advantages of working in the IMF have a 



transverse charge densities 



42J44 



44) and IPDs 



ready 



been emphasized in the derivation of 



13Ml6l] from form factors (FFs) and GPDs, 



respectively. Analogously, they will be exploited here to arrive at a definition of Wigner 
distributions which is not spoiled by relativistic corrections. 

Introducing two lightlike four- vectors n± satisfying n + • n_ = 1, we write the light- 
cone components of a generic four- vector a as [a + ,a~ ,a±] with a* 1 = a ■ n T . Similarly to 
30], we define the Wigner operators for quarks at a fixed light-cone time y + = 



Refs. 



29 



as follows 



W [r] (b ± ,k A 



x 



dz d 2 z± 



(27TJ 



V(y-f)r>W(y + 



i) 



)3 D ' ' r 1 2/| 2 +=0 

with = [0,0, b±], p + the average nucleon longitudinal momentum and x = k + /p + the 



average fraction of nucleoli longitudinal momentum carried by the active quark. The su- 
perscript T stands for any twist-two Dirac operator T = 7 + , 7 + 7s, ia^ +r j^ with j — 1,2. A 
Wilson line W = W(y — §, y + ||n) ensures the color gauge invariance of the Wigner oper- 



ator, connecting the points (y — | ) and 



and (y + f ) + oo • n by straight lines [31]. 



via the intermediary points — |) + oo • n 



We define the Wigner distributions in terms of the matrix elements of the Wigner oper- 
ators sandwiched between nucleon states with polarization S as follows 

p [T] (b±A,x,S) = J ^( P + ^,S\W^(b ± ,k ± ,x)\p\-^,S). (2) 

Thanks to the properties of the Galilean subgroup of transverse boosts in the IMF 15|,|45], we 
can form a localized nucleon state in the transverse direction, in the sense that its transverse 
center of momentum is at the position f± : 

\p + ,r±) = J ^ e " ¥± '" X b + '^)- ( 3 ) 

The Wigner distributions defined according to Eq. fl2]) can then be written in terms of these 
localized nucleon states as 

P [r] (b ± ,k ± ,x,S) = J d 2 D L (p + ,-^,S\W^(b ± ,k ± ,x)\p + ,fy,S), (4) 

where D± is the transverse distance between the initial and final centers of momentum. 
Note that the nucleon in our definition of the Wigner distributions has vanishing average 
transverse position and average transverse momentum, see Eqs. ()2]) and (j4j). This allows 
us to interpret the variables b± and k± as the relative average transverse position and the 
relative average transverse momentum of the quark, respectively. 
Using a transverse translation in Eq. (T5]), we find 

p^(b ± ,k ± ,x,S) = J ^e- iK ^wFKA ± ,k±,x,S), (5) 



where are the quark-quark correlators defining the GTMDs 3l| for A 



W^(A ± , fejL, x, S) = (p + , £s S\W^(0 ± , k ± , x)\p + , S) 

1 r J ~— A2 y _ _ _ (6) 

This means that the Wigner distributions defined as in Eq. ([2]) are the two-dimensional 
Fourier transforms of GTMDs, just like transverse densities and IPDs are two-dimensional 



Fourier transforms of FFs and GPDs, respectively. Contrarily to all the other distribu- 
tion functions, the GTMDs are in general complex-valued functions. However the two- 
dimensional Fourier transforms of the GTMDs are always real-valued functions, in accor- 
dance with their interpretation as phase-space distributions. 



B. Three-Dimensional Probability Densities 

Wigner distributions cannot have a strict probabilistic interpretation, since Heisenberg 
uncertainty relations prevent to determine at the same time position and momentum of a 
particle. Accordingly, Wigner distributions are not positive definite. Nevertheless, inte- 
grating out position and/or momentum variables, Wigner distributions reduce to genuine 
probability distributions. There are in particular four types of three-dimensional probability 
densities: 

• Integrating over b±_ amounts to set A i = i , and so the Wigner distributions reduce 



to the standard TMD correlators $^ 



31 



34] 



/ d 2 b ± p [r] (b ± , k ± , x, S) = W [r \6 ± , k ± , x, S) 
J (7) 

which can be interpreted as quark densities in three-dimensional momentum space; 
Integrating over k±_ amounts to set z± = 0j_, and so the Wigner distributions reduce 



to two-dimensional Fourier transforms of the standard GPD correlators 

12 



31 



34] 



/ 



with 



fF\K ± ,x,s) = IJ ^J xp+Z ~ <p + , ^,^l^(-f)rw^(f)b + ,-^,^>L += ^ =0 , 

(9) 

where the modulus of a general transverse vector a± is indicated as a±. In other words, 
one recovers the IPDs which can be interpreted as quark densities in the transverse 
position space and longitudinal momentum space; 
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• Integrating over b y and k x amounts to set = z x = 0, and so the Wigner distributions 
reduce to new three-dimensional quark densities 

j dbydk xP N(b ± ,k ± ,X,S) = p [r] (b x ,ky,X,S). (10) 

The variables b x and k y refer to two orthogonal directions in the transverse plane and 
so are not subjected to Heisenberg uncertainty relations; 

• Integrating over b x and k y amounts to set A x = z y = 0, and so the Wigner distributions 
reduce to other new three-dimensional quark densities 

j db X dky P [T] (b ± ,k ± ,X,S) = p [T] (by,k X ,xj). (11) 

There are a priori no simple relations between the quark densities in Eqs. (flUj) and 
( fTTT) . except when the quark and nucleon polarizations have no transverse components. 
In this case, the only privileged directions in the transverse plane are b± and k±, and 
we have p [r] (b ± , k±, x, e z ) = p [r] (-6±, k±, x, e z ) for V = 7 + ,7 + 7 5 . 



C. Wigner Distribution with Longitudinal Polarizations 



On the one hand, there are in total 16 GTMDs at twist-two level 31|. On the other 
hand, the quark and nucleon can be either unpolarized or polarized along three orthogonal 
directions, which means 16 configurations. All the 16 configurations can be written in terms 
of 16 independent linear combinations of the GTMDs. We will not present all of them in this 
study. To keep the discussion relatively simple, we focus on cases without any transverse 
polarization. 

The Wigner distribution of quarks with longitudinal polarization A in a nucleon with 
longitudinal polarization A is obtained for T = 7 + 1+ ^ 75 and S = Ae z 



pA\(b±,k ± ,x) = - p [7+] (6j_, k±,x,Ae z ) + Ap [7+7sl (6j_, k±,x,Ae z ) 



(12) 



We decompose it as follows 



PA\(b_i_, k L ,x 
1 

~ 2 



Puu{b±,k±,x) + Ap LU (b ± ,k ± ,x) + \p UL (b±,k ± ,x) + A\p LL (b ± ,k±,x) , (13) 
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where 



Puu(b±,kx,x) = ^ P b+] (b ± ,k±,x,+e z ) + p [7+1 (6j_, -e 2 ) (14) 
is the Wigner distribution of unpolarized quarks in an unpolarized nucleon; 

PLu(b.L,k±,x) = - p lj+] (b ± , k±,x, +e z ) - p b+] (b ± , k±,x, -e z ) (15) 

represents the distortion due to the longitudinal polarization of the nucleon; 

p UL (bx, k±, x)= l - J> + ^ (6 ± , k ± , x, +e z ) + p^ (b L , k ± , x, -e z )] (16) 
represents the distortion due to the longitudinal polarization of the quarks, and 

p LL (b±,k ± ,x) = ~ [p [ ^ ] (b ± ,k ± ,x,+e z ] -p [ ^' 5] (b ± ,k ± ,x,-e z )] (17) 

represents the distortion due to the correlation between quark and nucleon longitudinal 
polarizations. These four contributions can be written as 

Puu(b±, k±, x) = 0, k\, k ± ■ b ± , b] 

PLu(b±,k±,x) -- 



I d - - - -. 

e?fc l i — r Fi^ix, 0, kj_, k± ■ b±, b 2 ± 



M 2 



db\ 



- -> 1 ■ ■ d - - - 

Pudb±, k ± , x) = — € l k ±ggr Gi t i{x, 0, k 2 x , k x ■ b±, b 

PLL(b±, k ± , x) = GiAx, 0, k 2 ± , k ± ■ b±, b^), 



I), 



(18a) 
(18b) 

(18c) 
(18d) 



where the distributions X = T\^T\^ G\,\, £1,4 are the Fourier transforms of the correspond- 
ing GTMDs X = F 1A , F 1A , G 1}1 , G 1A introduced in Ref. |31| 

d 2 A 



X(x,£,k ± ,k± ■ b ± ,bj 



-iA 1 b 



(27T) 5 



X(x,£,ki,k ± - A ± ,A 



(19) 



In Eq. (JT8"]) the two-dimensional antisymmetric tensor has been used with e 
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,21 



M is the nucleon mass and roman indices are to be summed over. Integrating out b± 
or k± kills the contributions plu and pul, showing that there exists no TMD or GPD 
corresponding to -F 14 and Gi A . These GTMDs carry therefore completely new information 
about the nucleon structure. On the other hand, the contributions puu an d Pll survive 
both integrations. It follows that the GTMD F11 can be seen as the mother distribution of 
the TMD fx and the GPD H 

A(x, k\) = [ d 2 b ± F 1A (x, 0, k 2 ± , k ± ■ 6 X , bl) = F 1A (x, 0, k 2 ± , 0, 0), 



H(x,0,A 2 



-5 ^±Ji 



(20a) 
(20b) 



and the GTMD G\ t 4 as the mother distribution of the TMD gn and the GPD H 



g 1L (x, %l) = J dV Gi,a(x, 0, k 2 ± , k ± ■ b ± , b\) = G 1>4 (x, 0, k\, 0, 0), (20c) 

H(x, 0, A 2 ± ) = J d 2 k ± G 1A (x, 0, k 2 ± , k ± ■ A ± , A 2 ± ). (20d) 
Integrating out all the variables, one naturally gets 

J dx d 2 k ± d 2 b ± pjj^, k ± , x) = N q , (21a) 

J dxd 2 k ± d 2 b ± p g LU (b ± ,k L ,x) = 0, (21b) 

J dxd 2 k ± d 2 b ± p q UL (b ± ,k ± ,x) = 0, (21c) 

/ dxd 2 k ± d 2 b ± p q LL (b ± ,k ± ,x) = Aq, (21d) 



where the index q indicates the contribution of the quark of flavor g, N q is the valence-quark 
number (N u = 2 and N d = 1 in the proton) and Aq is the axial charge. Note that Eq. (121 bp 
tells us that the valence-quark number does not depend on the nucleon polarization and 
Eq. (I21cl) means that in an unpolarized nucleon there is no net quark polarization. 



D. Quark Orbital Angular Momentum 

Quantifying quark orbital angular momentum (QAM) inside the nucleon is essential in 
order to solve the so-called "spin crisis", see e.g. Q]- Almost 15 years ago, Ji derived a 
sum rule that allows one to extract the total quark contribution to the nucleon spin from a 



combination of GPDs 48| 



J* = ^Jdxx [H q (x, 0, 0) + E q (x, 0, 0)] . (22) 

By subtracting half of the axial charge Aq = j dxH q (x,0,0) which is interpreted as the 
spin contribution of quarks with flavor q to the nucleon spin, one gets the quark OAM 
contribution 

L q z = ^Jdx{x [H q (x, 0, 0) + E q (x, 0, 0)] - H q {x, 0,0)}. (23) 

From a density point of view, this result is surprising in the sense that the extraction of 
the quark OAM along the z-axis involves the GPD E which appears only in a transversely 
polarized nucleon. Note however that E describes the amplitude where the nucleon spin 



flips while the quark light-cone helicities remain unaffected, implying therefore a change by 
one unit of OAM between the initial and final nucleon states. 

More recently it has been suggested, based on some quark models, that the TMD h^ T 



may also be related to the quark OAM 



49-152 1 



C\ = - I dxd 2 k ± ^- 2 hf^x,kl). (24) 



Note that one expects in general C q z ^ L\ in a gauge theory, see e.g. [53(. Once again, 
from a density point of view, this expression is surprising in the sense that it involves the 
TMD hjj, which describes the distribution of transversely polarized quarks in a transversely 
polarized nucleon. Note however that h^ T corresponds to the amplitude where the nucleon 
and active quark longitudinal polarizations flip in opposite directions, involving therefore a 
change by two units of OAM between the initial and final nucleon states. 

Clearly, Wigner distributions provide much more information than GPDs and TMDs 
as they contain also the full correlations between quark transverse position and three- 
momentum. Furthermore, once the Wigner distributions are known, it is rather straightfor- 
ward to compute physical observables. One has just to take the phase-space average as if 
the Wigner distributions were classical distributions 

(A)V\S) = J dxd 2 k ± d 2 b ± A(b ± ,k ± ,x)p [r] (b ± ,k ± ,x,S). (25) 

In particular, we can write the average quark OAM in a nucleon polarized in the ^-direction 

as 

i\ = (L«) [7+] (e,) = J dxd 2 k ± d 2 b ± (b ± x kA_) z P [ i +]q (b x ,k ± ,x,e z ) 

= J dxd 2 k ± d 2 b ± (b ± x [p q uu (b±,k ± ,x) + p q LU (b ± ,k ± ,x)]. (26) 

From Eq. ( 118al) . it is clear that 

J dxd 2 k ± d 2 b ± (b ± x k ± ^ p^(6±,fc±,x) = 0, (27) 
which means that in an unpolarized nucleon there is no net quark OAM 1 . Using now 



1 An unpolarized nucleon has no spin, which means that the total quark and gluon angular momentum 
contributions have to sum up to zero. By rotational invariance, one expects all the four contributions 
(spin and OAM of quarks and gluons) to vanish identically. The angular momentum sum rule for an 
unpolarized nucleon is therefore trivially satisfied. 
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Eq. (118bj) and integrating by parts, we find that the quark OAM £ q z reads 



(28) 



An interesting issue which deserves further investigation is the relation between L q z in 
Eq. fl23|) and £ q z in Eq. ( 1281) . As discussed in the following sections, in models without 
gauge-field degrees of freedom one finds that the two definitions give the same results for 
the total quark contribution to the OAM, but not for the separate quark-flavor contributions. 
However, this remains to be confirmed in more complex systems, when the contribution of 
the Wilson line is explicitly taken into account. 

Wigner distributions allow us also to study the correlation between quark spin and OAM, 
which we define as 



CI 



dxd k± d b± [b± x kj_j Pi; L (b±, k±, x) 
dxd 2 k x ^G q 1A (x, 0,k 2 ± , 0,0), 



(29) 



where we have used Eq. f)18cl) . For CJ > the quark spin and OAM tend to be aligned, 
while for C\ < they tend to be antialigned. Finally, note that from Eq. f!18dj) one has 



dx d 2 k± d 2 b± [b± x kj_j PLL(b±,k±, x) = 0, 
reflecting like Eqs. f!21bj) . f)21cl) and fl27|) the isotropy of space. 



(30) 



III. RESULTS AND DISCUSSIONS 



Since it is not known how to extract Wigner distributions or GTMDs from experiments, 
one has to rely on phenomenological models. We studied the Wigner distributions in the 
light-cone constituent quark model (LCCQM) 32|, |33|, |35| and the light-cone version of the 



chiral quark-soliton model (xQSM) restricted to the three-quark sector [38|, |40|, |4l|, |54], l55| . 
using the general formalism developed in Ref. 34j for the overlap representation of the quark- 



quark correlator in terms of light-cone wave functions. We neglect the contribution from 
gauge degrees of freedom, and in particular from the Wilson line in the Wigner operator ([1]). 
As the resulting distributions are very similar in both models, we will present only those 
from the LCCQM. However, when discussing more quantitative aspects, we will also report 
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the numerical values from the xQSM. Furthermore, we will discuss only the first x moment 
of the Wigner distributions 

p{b±,k±) = J dxp(b ± , k±,x), (31) 

i.e. purely transverse four- dimensional phase-space distributions (two transverse position 
and two transverse momentum coordinates), referred to as transverse Wigner distributions. 



A. Unpolarized Quarks in an Unpolarized Nucleon 

We start the discussions with puu(b±, k±), the transverse Wigner distribution of unpo- 
larized quarks in an unpolarized proton. In Fig. [T] we present the distributions in impact- 
parameter space with fixed transverse momentum k± = k± e y and k± = 0.3 GeV (upper 
panels), and compare them with the distribution in transverse- momentum space with fixed 
impact parameter b± = bj_e y and b± = 0.4 fm (lower panels). The left (right) panels re- 
fer to the u (d) quarks. We observe a distortion in all these distributions which indicates 
that the configuration bj_ _L kj_ is favored with respect to the configuration bj_ \\ kj_. This 
can be understood with naive semi-classical arguments. The radial momentum (k± ■ b±) b± 
(b± = b±/b±) of a quark is expected to decrease rapidly in the periphery because of confine- 
ment. The polar momentum k± — (k± ■ b±) b± receives a contribution from the orbital motion 
of the quark which can still be significant in the periphery (in an orbital motion, one does 
not need to reduce the momentum to avoid a quark escape). This naive picture also suggests 
that this phenomenon should be stronger as we go to peripheral regions (b± and to high 
quark momenta (k± ^>). Such a behavior is indeed observed in our model calculations and 
can be quantified in terms of the average quadrupole distortions Q l ^(k±) and Q^(b±) defined 
as 

" V > fd% x bl Puu (b ± ,k x ) 

.. ^ x , s / .... A f d 2 k L (lk\k\ -S^k 2 ,) PuU (b ± ,k L ) 

QtibJ = Q k (b x ) [2VV - S») = J - ^-rr- 2 f " \ (32b) 

v J J d 2 k ± kj_p uu {b ± ,k ± ) 

where i,j = x,y. The distortions calculated in the LCCQM are tabulated in Tabled! We 
note that the quadrupole distortions of u and d quarks are very similar. For increasing 
values of k± and b± the distortions get more and more pronounced, and at the same time 
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Pvv W(GeV*-fm>)] p d vv [l/(GeV*-fm>)] 




-0.4 -0.2 0.0 0.2 0.4 -0.4 -0.2 0.0 0.2 0.4 

k x [GeV] k x [GeV] 



FIG. 1: The transverse Wigner distributions of unpolarized quarks in an unpolarized proton. Upper panels: 
distributions in impact-parameter space with fixed transverse momentum k± = k±_ e y and k± = 0.3 GeV. 
Lower panels: distributions in transverse-momentum space with fixed impact parameter b± = b±^ e y and 
b± = 0.4 fm. The left (right) panels show the results for u (d) quarks. 

the spread of the distributions shrinks towards the center. From Fig. [T] we also note that 
the spread of the distributions is smaller for u quarks than for d quarks, especially in the 
transverse-coordinate space. This reflects the fact that u quarks are more concentrated at 
the center of the proton, while the <i-quark distribution has a tail which extends further at 
the periphery of the proton. 

From Eq. ( jl8a|) we see that puu(b±, k±) = puu(b±, k±, k±-b±). This explains the left-right 
symmetry in Fig. [T] and implies that the quark is as likely to rotate clockwise as to rotate 
anticlockwise. In Fig. [1] we also observe a top-bottom symmetry. Such a symmetry is not 
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TABLE I: The average quadrupolc distortions of the transverse Wigner distribution of unpolarizcd quarks 
in an unpolarized proton from the LCCQM. Sec Eq. for the definition of Qb(k±) and Qk(b±)- 



k x [GeV] 


Qb( 

u 


k±) 
d 


b± [fm] 


Qk(b±) 

n, \ _l_ / 

u d 

















0.1 


-0.04 


-0.03 


0.2 


-0.15 -0.10 


0.2 


-0.14 


-0.13 


0.4 


-0.24 -0.19 


0.3 


-0.30 


-0.27 


0.6 


-0.29 -0.27 



a general property of the Wigner distribution puu-, but follows from the fact that in our 
calculations there are no explicit gluons. Indeed, time-reversal invariance implies that the 
real part of the GTMDs is T-even (e) while the imaginary part is T-odd (o). Hermiticity 
tells us that the four GTMDs X = F 14 , G^i, Gi j4 satisfy the relations 

X e (x, £, hi, k± ■ A x , Ai) = X e (x, k 2 ± , -k ± ■ A ± , A 2 ± ), 
X°(x, £, P ± , k ± ■ A ± , Ai) = -X°(x, P ± , -k ± ■ A ± , Ai). 

This means that for £ = 0, X e is an even function of k± ■ A± while X° is an odd function of 
k± ■ Aj_. It follows that the Fourier transforms X of these GTMDs with respect to A^ are 
real-valued functions. The contribution X e is an even function of k± ■ b± while X° is an odd 
function of k^-b±. Since we have no explicit gluons and therefore no final-state interactions, 
our GTMDs are real. It follows that X = X e explaining the top-bottom symmetry of Fig. [TJ 
The dominant effect of final-state interactions would be to shift up or down the distributions. 

As mentioned earlier, Wigner distributions have only a quasi-probabilistic interpretation 
due to Heisenberg uncertainty relations. A genuine probabilistic interpretation can be re- 
covered only when integrating out certain variables. If we integrate out b± or k±, we reduce 
to the unpolarized TMD and GPD, respectively. In these cases, the distortion we observed 
in the Wigner distribution puu is completely absent and we are left with axially symmetric 
distributions, see Eq. ( 1201) . By integrating over one momentum and one coordinate variables 
which are not conjugated, we obtain the probability densities p of Eqs. ( flOl) and ( flTT) . In 
Fig. |2]we show the probability density puu(b x ,k y ) which gives the correlation between b x 
and k y . We observe that puu(b x , k y ) is maximum at the center, b x = k y = 0, and decreases 
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FIG. 2: The mixed transverse densities p{b x ,k y ) of unpolarized u quarks (left panel) and unpolarized d 
quarks (right panel) in an unpolarized proton. 

in the outer regions of the phase space, with the equidensity lines in each quadrant of Fig. [2] 
having approximately a linear dependence in b x and k y . Furthermore, we clearly see that 
the width of the densities in k y is similar for u and d quarks, while it is more extended in b x 
for d quarks than for u quarks. 



B. Unpolarized Quarks in a Longitudinally Polarized Nucleon 

We now consider piu(b±, k±), the distortion of the transverse Wigner distribution of un- 
polarized quarks due to the longitudinal polarization of the proton. In Fig. [31 the upper 
panels show the distortions in impact-parameter space for u (left panels) and d (right panels) 
quarks with fixed transverse momentum k± = k± e y and k±_ = 0.3 GeV, while the lower pan- 
els give the corresponding distortions in the transverse-momentum space with fixed impact 
parameter b± = b±e y and b± = 0.4 fm. We observe a clear dipole structure in both these 
distributions, with opposite sign for u and d quarks. The corresponding distortions of the 
mixed transverse densities p(b x ,k y ) are shown in Fig. H]for the u (left panel) and d (right 
panel) quarks. In this case, we observe a quadrupole structure. These multipole structures 
are due to the explicit factor eY^l^§~ m -^q- f!18bj) which breaks the left-right symmetry 
in Fig. [3] allowing therefore non- vanishing net OAM. We learn from these figures that the 
OAM of u quarks tends be aligned with the nucleon spin, while the OAM of d quarks tends 
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FIG. 3: The distortions of the transverse Wigner distributions of unpolarized quarks due to the spin of 
the proton (pointing out of the plane). Upper panels: distortions in impact-parameter space with fixed 
transverse momentum k±_ = k± e y and k± = 0.3 GeV. Lower panels: distortions in transverse-momentum 
space with fixed impact parameter b± = b± e y and b± = 0.4 fm. The left (right) panels show the results for 
u (d) quarks. 

be antialigned with the nucleoli spin. In particular, we notice that the distortion induced 
by the quark OAM is stronger in the central region of the phase space (A;^ <^ and b± <C), 
for both u and d quarks. The distortion in the b± space (see upper panels of Fig. [3]) is 
more extended for d quarks than for u quarks, whereas the opposite behavior is found for 
the distortion in the k± space (see lower panels of Fig. [3]). In the case of d quarks, we also 
observe a sign change of the distributions in the outer regions of phase space (k± 3> and 
b± 3>) which corresponds to a flip of the local net quark OAM. 

According to Eq. (|28|) . the integral over the phase space of the distribution in Fig. [3] 
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FIG. 4: The distortions of the mixed transverse densities p(b x , k y ) of unpolarizcd u quarks (left panel) and 
d quarks (right panel) due to the spin of the proton (pointing out of the plane) . 



TABLE II: The results for quark orbital angular momentum (see Eqs. (|23|) . ([24]) and (|28[0 and anomalous 
magnetic moment n q from the LCCQM and the xQSM for u-, d- and total (u + d) quark contributions. 





Model 




LCCQM 






XQSM 






Q 


u 


d 


Total 


u 


d 


Total 


a 


Eq. d28]) 


0.131 


-0.005 


0.126 


0.073 


-0.004 


0.069 


Li 


Eq. d23D 


0.071 


0.055 


0.126 


-0.008 


0.077 


0.069 


a 


Eq. {MD 


0.169 


-0.042 


0.126 


0.093 


-0.023 


0.069 




1.867 


-1.579 


0.288 


1.766 


-1.551 


0.215 



multiplied by (b± x k±) z gives the expectation value of the quark OAM. The corresponding 
results for u-, d- and total (u + d) quark contributions are reported in the first row of Table ILT1 
We give the predictions from both the LCCQM and the xQSM. They are compared with 
the corresponding results for the orbital angular momentum L q z and C q z obtained from the 
definitions in Eqs. and (T24"|) . respectively. We note that all the three definitions give the 
same results for the total quark OAM, but differ for the separate quark-flavor contributions. 
In agreement with the interpretation of the TMD results given in Ref . 34J , there is more net 
quark OAM in the LCCQM (J2 q L q z = 0.126) than in the xQSM (J2 q L l = 0.069). For the 
individual quark contributions, both the LCCQM and the xQSM predict that l\ and C q are 
positive for u quarks and negative for d quarks, with the u-quark contribution larger than 
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the ci-quark contribution in absolute value. For L\ the LCCQM predicts the same positive 
sign for the u and d contributions, with the isovector combination L u z — L d > 0, similarly 
to a variety of relativistic quark model calculations 53] . Instead, the yQSM gives L" < 
and L d > 0, and therefore L u z —L d z < 0, in agreement with lattice calculations 0, Note 
however that the quark angular momenta J| of Eq. (122j) are very similar in both models. 
We find J u z = 0.569 (0.566) and J d z = -0.069 (-0.066) in the LCCQM (xQSM). This is due 
to the fact that the difference in L q z between the two model predictions is compensated by 
the different results for the spin contribution, namely Am = 0.995 (1.148) and Ad = —0.249 
(-0.287) in the LCCQM (xQSM). 

In the last row of Table HT1 we also give the results for the quark anomalous magnetic mo- 
ment K q which is intimately connected to quark OAM. In particular, it is well known within 
the light-cone wave function description of hadrons that a state can have anomalous mag- 



netic moment only in the presence of nonzero OAM components 58j-[6l|. Furthermore, the 



anomalous magnetic moment gives a measurement of the correlation between the transverse 
spin of the nucleon and the orbital motion of quarks, as observed in the IPDs for unpo- 

nn 

larized quarks in a transversely polarized nucleon [14J, [15[. We find the following pattern 
< —n d < k u which coincides with the ones for C q , and if but not for L%. 



C. Longitudinally Polarized Quarks in an Unpolarized Nucleon 

We now discuss puL(b±, fcjj, the distortion of the transverse Wigner distribution due to 
the longitudinal polarization of quarks in an unpolarized proton. In Fig. [5] we show the 
distortions, both in impact-parameter space with fixed transverse momentum kj_ = k± e y 
and k± = 0.3 GeV (upper panels) and in transverse-momentum space with fixed impact 
parameter b± = b± e y and b± = 0.4 fm (lower panels), of the distributions of u (left panels) 
and d quarks (right panels). The corresponding distortions of the mixed transverse densities 
p(b x , k y ) are shown in Fig. [6], for u (left panel) and for d (right panel) quarks. 

Like in the case of the pm distributions, the dipole and quadrupole structures are due to 
the explicit factor ej^l^fr in Eq. f|18c[) . We learn from these figures that the quark OAM 
and the quark spin tend to be aligned for both u and d quarks. The size of the distributions 
are similar for u and d quarks. However, since there are effectively twice more u quarks than 
d quarks in the proton, the alignment is more pronounced for the d quarks than the u quarks. 
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FIG. 5: The distortions of the transverse Wigner distributions due to the spin of the quarks (pointing out 
of the plane) in an unpolarized proton. Upper panels: distortions in impact-parameter space with fixed 
transverse momentum k±_ = k± e y and k± = 0.3 GeV. Lower panels: distortions in transverse-momentum 
space with fixed impact parameter b± = b± e y and b± = 0.4 fm. The left (right) panels show the results for 
u (d) quarks. 

The correlation CJ between quark spin and OAM in the e z direction can be calculated using 
the definition in Eq. ( 1291) . The results for both the LCCQM and the xQSM are given in 
Table Hill As anticipated, we find > and C d z > 0, with larger values in the LCCQM 
than the xQSM. 

In Table IIHI we also give the results for the tensor anomalous magnetic moment k\ which 
measures the correlation between the transverse spin and the transverse OAM of the quark 
in an unpolarized nucleon, as observed in the IPDs for transversely polarized quarks in an 
unpolarized nucleon 13, |62|. We find > > 0, which coincides with the pattern of C\. 
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FIG. 6: The distortions of the mixed transverse densities p(b x ,k y ) due to the spin (pointing out of the 
plane) of u quarks (left panel) and d quarks (right panel) in an unpolarizcd proton. 

TABLE III: The results for quark spin-OAM correlation Cf (see Eq. (12911 ) and anomalous tensor magnetic 
moment nip obtained in the LCCQM and the xQSM for u-. d- and total (u + d) quark contributions. 



Model 

q 


LCCQM 
u d Total 


XQSM 
u d Total 


CI Eq. ([29D 


0.227 0.187 0.414 


0.130 0.109 0.239 


Kirp 


3.947 2.581 6.528 


3.832 2.582 6.414 



However, at variance with Cf , the values for k\ from the LCCQM 33( and the xQSM 
are very similar. 



D. Longitudinally Polarized Quarks in a Longitudinally Polarized Nucleon 

We proceed with the discussion of PLi{b±, k±), the distortion of the transverse Wigner 
distribution due to the correlation between the longitudinal polarizations of the quarks and 
the proton. In Fig. [7] we show the distortions, both in impact-parameter space with fixed 
transverse momentum k± = k± e y and k± = 0.3 GeV (upper panels) and in transverse- 
momentum space with fixed impact parameter b± = b± e y and b± = 0.4 fm (lower panels), of 
the distributions of u (left panels) and d quarks (right panels). The corresponding distortions 
of the mixed transverse densities p(b x , k y ) are shown in Fig. [HJ As one already knows from 
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FIG. 7: The distortions of the transverse Wigner distributions due to the correlation between the quark 
spins and the proton spin (pointing out of the plane). Upper panels: distortions in impact-parameter space 
with fixed transverse momentum k± = k±_ e y and fcj_ = 0.3 GeV. Lower panels: distortions in transverse- 
momentum space with fixed impact parameter b± = b± e y and b± = 0.4 fm. The left (right) panels show 
the results for u (d) quarks. 

the axial charges, the u-quark polarization tends to be parallel to the nucleon spin, while the 
(i-quark polarization tends to be antiparallel. Accordingly, the distributions are positive for 
u quarks and negative for d quarks. The new information is about the distribution in phase 
space of these polarizations (see Eq. (121d[) ) . It appears that the quark polarization receives 
its main contribution from the central region of the phase space. Interestingly, the average 
quark polarization changes sign for sufficiently large bj_ or k±, preferably when bj_ and k± 
are aligned (see Fig. [7]). From Eq. ( 118dl) we see that Pll(6j_, k±) = PLi(b±, k±, k± ■ bj_), 
explaining the left-right symmetry in Fig [71 It follows that pn cannot contribute to the net 
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FIG. 8: The distortions of the mixed transverse densities p(b x , k y ) due to the correlation between the spin 
of u quarks (left panel) and d quarks (right panel) , and the proton spin (pointing out of the plane) . 

quark OAM, as required by the isotropy of space (see Eq. (1501) ). 

Combining as in Eq. f p~3|) the Wigner distribution of unpolarized quarks in an unpolarized 
proton puu with the distortions plu, Pul and pll, we obtain the Wigner distribution p\\ of 
longitudinally polarized quarks in a longitudinally polarized proton. In Fig. [9J the transverse 
Wigner distributions of u and d quarks with polarization A =t, 4 in a proton with polar- 
ization A =f are shown in the impact-parameter space with fixed transverse momentum 
k± = k±e y and k± = 0.3 GeV. The corresponding mixed transverse densities p~A\(b x ,k y ) 
are shown in Fig. [TpJ The deformation induced by the quark and nucleon polarizations 
is clearly visible in the sideway shifts of the distributions in Fig [9j In particular, when 
the quark and nucleon polarizations are parallel (antiparallel) the shift is in the positive 
(negative) b x direction, see upper (lower) panels. 

We learned from puj and pul that the u-quark OAM tends to be aligned with both the 
quark and proton polarizations. When the u quark has polarization parallel to the nucleon 
spin, the contributions pm and pm interfere constructively resulting in a sideway shift in the 
positive b x direction. When the u quark has polarization antiparallel to the nucleon spin, the 
contributions piu an d pm interfere destructively. Since the correlation between the OAM 
and the quark spin is stronger than the correlation between the OAM and the nucleon spin 
(see Figs. [3]and|5]), it results a sideway shift in the negative b x direction. For the d quark, 
we learned from p LU and pul that the OAM tends to be aligned with the quark polarization 
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FIG. 9: The transverse Wigner distributions of longitudinally polarized quarks in a longitudinally polarized 
proton (A =f pointing out of the plane) in impact-parameter space with fixed transverse momentum k±_ = 
fc_L e y and k± = 0.3 GeV. Upper panels: distributions of quarks with polarization parallel to the nucleon 
spin (A =t). Lower panels: distributions of quarks with polarization antiparallel to the nucleon spin (A —l). 
The left (right) panels show the results for u (d) quarks. 

but antialigned with the nucleon polarization. When the d quark has polarization parallel to 
the nucleon spin, the contributions piu and pui interfere destructively. Once again, we can 
see from Figs. [3] and [5] that the correlation between the OAM and the quark spin is stronger 
than the correlation between the OAM and the nucleon spin, resulting in a sideway shift 
in the positive b x direction. When the d quark has polarization antiparallel to the nucleon 
spin, the contributions piu and pui interfere constructively resulting in a sideway shift in 
the negative h x direction. The sideway shifts are more apparent in p"^ and because 
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FIG. 10: The mixed transverse densities pA\{bx,k y ) of longitudinally polarized quarks in a longitudinally 
polarized proton (A =j" pointing out of the plane). Upper panels: probability densities of quarks with polar- 
ization parallel to the nucleon spin (A =t)- Lower panels: probability densities of quarks with polarization 
antiparallel to the nucleon spin (A =\). The left (right) panels show the results for u (d) quarks. 

the contributions puu and Pll partially cancel in these cases. In an analogous way we can 
understand the densities in the (b x , k y ) plane shown in Fig. [TOJ 



IV. CONCLUSIONS 



In this work we presented a study of the quark Wigner functions which provide the 
full phase-space description of the quark distributions in the nucleon. Using the light- 
front formalism, we derived the Wigner distributions as two-dimensional Fourier transforms 
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of the GTMDs from the transverse- momentum transfer Aj_ to the impact parameter b±. 
Therefore these distributions provide us with images of the nucleon in five dimensions, 
namely two position and three momentum coordinates. This derivation is not spoiled by 
relativistic corrections and is completely analogous to the interpretation of transverse charge 
densities and impact-parameter dependent parton distributions as two-dimensional Fourier 
transforms of form factors and GPDs, respectively. However, Wigner distributions can not 
have a strict probabilistic interpretation, since Heisenberg uncertainty relations forbid to 
localize a particle and to determine its momentum at the same time. Accordingly, the 
Wigner distributions are not positive definite. Only in particular limits one can recover a 
density interpretation. This is the case for the known projections of the Wigner distributions 
to the three-dimensional densities in the momentum or in the impact-parameter space at 
fixed x, corresponding to TMDs and IPDs, respectively. On top of them, we introduced two 
new types of densities mapping the nucleon as functions of one spatial and one momentum 
variables in the transverse plane which are not conjugated, and therefore not constrained 
by the uncertainty principle. 

In general the GTMDs are complex-valued functions. However their two-dimensional 
Fourier transforms are always real-valued functions, in accordance with the interpretation of 
the Wigner distributions as phase-space distributions. At leading twist, there are 16 Wigner 
distributions which can be disentangled by varying the nucleon and quark polarizations. 
We focused on the four cases without transverse polarizations, namely the distributions 
of unpolarized/longitudinally polarized quarks in an unpolarized/longitudinally polarized 
nucleon. Furthermore, we considered only the quark contribution, neglecting all the gauge- 
field degrees of freedom. In this case, the imaginary part of the GTMDs is zero, and by 
hermiticity and time-reversal invariance the corresponding Wigner distributions are even 
functions of k± ■ b±. 

The results for the Wigner distributions within a light-cone constituent quark model and 
the light-cone version of the chiral quark-soliton model are very similar, and allowed us to 
sketch some general features about the behavior of the quarks in the nucleon when observed 
in the b± plane at fixed k±, or in the k± plane at fixed b±. In particular, the Wigner 
distributions of unpolarized quarks in an unpolarized nucleon are not axially symmetric. 
This deformation can be explained with naive semi-classical arguments, as a consequence of 
confinement which limits the radial motion of the quark with respect to its orbital motion. 
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The effect becomes more pronounced at larger values of k± (b±) in the impact-parameter 
(momentum) space, as found by calculating the corresponding quadrupole distortions. Fur- 
thermore, we observed that the spread of the distributions is smaller for u quarks than for 
d quarks, especially in the b± space, revealing that the u quarks are more concentrated at 
the center of the proton, while the d-quark distribution has a tail which extends further at 
the periphery of the proton. 

The case of unpolarized quarks in a longitudinally polarized nucleon is particularly in- 
teresting because it allows us to calculate the phase-space average of the quark OAM. The 
corresponding values within the LCCQM and the light-cone xQSM have been compared 
with the results for the OAM from the Ji's sum rule, and from the definition in terms of 
the hyp TMD. We found that in models without gauge fields these three definitions give the 
same values for the total quark contribution to the OAM, while they differ for the individual 
u and d quark-flavor contributions. A peculiar result of the light-cone xQSM is that the 
isovector combination (u — d) of the OAM calculated from the Ji's sum rule is found to be 
negative, in agreement with lattice calculations and at variance with most of quark models. 

The distortion due to the longitudinal polarization of quarks in an unpolarized nucleon 
allowed us to study the correlation between the quark spin and OAM. This correlation, 
taking into account the effective number of u and d quarks, has been found to be stronger 
for d quarks than for u quarks. The same behavior has also been observed for the values of the 
tensor anomalous magnetic moments which measure the correlation between the transverse 
spin and the transverse OAM of quarks in an unpolarized nucleon. 

In the case of the distortion due to the correlation between the quark and nucleon spins, 
we were able to study the distributions of the axial charge in the phase space. They are 
positive for u quarks and negative for d quarks, and receive the main contribution from 
the central region of the phase space. An interesting finding is the fact that the average 
polarization of d quarks changes sign for sufficiently large b± and k±, preferably when b± 
and k± are aligned. 

Finally, taking into account all the four contributions discussed above, we visualized 
the combined effects induced on the distributions by the longitudinal polarizations of the 
quarks and nucleon. In all the examples we have studied, the prominent role of OAM and 
its correlation with quark and nucleon polarizations in shaping the quark distributions in 
the phase space has clearly emerged. Besides specific features related to the quark models 
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we used for discussing the results, we tried to emphasize the physical content of the Wigner 
distributions, and in particular the new information encoded in these distributions. Further 
studies of the Wigner distributions in different theoretical models can provide new insights 
on the quark and gluon dynamics. On the other hand, although Wigner distributions are not 
directly measurable from experiments, any new information on GPDs and TMDs extracted 
from experiments can be used to further constrain the multidimensional image of the nucleon 
in the quantum phase space. 
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